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SECTION  I 


INTRODUCTION 


The  "’•ohanical  properties  of  heavy  metal  alloys  avail- 
able for  use  as  long  rod  kinetic  energy  penetrators  vary  over 
a wide  range  of  values  depending  upon  the  constituents  of  the 
alloy  and  its  heat  treatment.  The  penetrator  designer  con- 
fronted with  such  choices  of  material  properties  must  attampt 
to  choose  the  particular  alloy  that  will  optimize  penetrability. 
To  accomplish  this,  he  must  understand  the  penetration  process 
and  hew  it  is  affected  by  penetrator  configurational  and  mate- 
rial parameters. 

Recently  two  theoretical  investigations  were  undertaken 
for  the  purpose  of  aiding  the  designer  in  his  choice  of  mate- 
rial properties  for  long  rod  penetration.  In  the  first  inves- 
tigation (Reference  1)  several  heavy  metal  penetrators  of 
different  designs  were  ranked  according  to  their  penetrability 
when  impacted  against  a 5 ,08-cm-thick  steel  target,  in  the 
second  effort  (Reference  2)  one  of  the  heavy  metal  penetrator 
designs  was  chosen  and  the  effect  of  varying  penetrator  mate- 
rial properties  was  examined.  In  particular , calculations 
involving  the  variation  in  penetrator  material  yield  strength. 


L-J.V-i  T m-rT* 
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and  the  effects  of  varying  these  properties  on  the  projectile's 
ability  to  penetrate  were  investigated.  The  second  effort  led 
to  several  tentative  conclusions  regarding  the  dependency  of 
final  crater  geometry  on  the  various  penetrator  material 
properties . 

The  current  effort  was  undertaken  in  order  to  substan- 
tiate or  otherwise  modify  some  of  the  earlier  conclusions  as 
well  as  to  investigate  ranges  of  parameters  not  considered  in 
the  previous  investigations.  In  this  report,  data  from  the 
current  effort  are  pooled  with  those  generated  in  the  previous 
investigations  and  the  analyses  reported  here  cover  all  of  the 
data  from  all  three  of  the  efforts. 

Section  II  of  this  report  briefly  describes  the  numer- 
ical and  material  models  employed  in  the  HELP  calculations. 


1.  Sedgwick,  R.  T.  and  L.  J.  Walsh,  "High  Density  Penetrator 
Performance  in  a Hard  Target,"  AFATL-TR-74-201,  ADB004857L, 
Air  Force  Armament  Laboratory,  Eglin  Air  Force  Baso, 
Florida,  December  1974. 

2.  Hageman,  L.  J.  and  R.  T.  Sedgwick,  "A  Parametric  Investi- 
gation of  the  Effects  of  Varying  the  Material  Properties 
of  a Heavy  Metal  Penetrator,"  AFATL-TR-75-152,  Air  Force 
Armament  Laboratory,  Eglin  Air  Force  Base,  Florida, 
November  1975. 
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In  Section  III  the  perimetric  studies  are  summarised  and  the 
calculation  matrix  presented,  section  IV  diacuesee  changes 
in  the  numerical  model  and  validation  of  these  changes  by 
comparison  with  experimental  data.  In  Section  V the  results 
from  the  calculations  are  presented  and  discussed.  These 
results  include  plots  which  show  the  dependency  of  final 
crater  depth  and  radius,  final  projectile  length,  and  energy 
coupled  to  the  target  on  such  material  properties  as  shear 
yield  strength,  failure  threshold,  shear  modulus. and  bulk 
modulus.  Section  Vl  presents  the  conclusions  drawn  from  the 
combined  current  and  previous  investigations  and  provides 
recommendations  for  additional  work. 

The  results  from  the  combined  parametric  investiga- 
tions of  the  effects  of  the  material  properties  of  heavy 
metal  penetrators  on  the  dimensions  of  the  crater  in  a steel 
target  indicate  that  available  numerical  techniques  can  pro- 
vide invaluable  guidance  to  the  penetrator  designer  regarding 
the  choice  of  penetrator  material  properties  for  a given 
penetrator  design  and  given  impact  conditions. 
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SECTION  II 


i 

i 


I 

NUMERICAL  AND  MATERIAL  MODELS  I 


2.1  NUMERICAL  TECHNIQUES 

In  the  present  material  parametric  study  of  heavy  metal 
penetrators,  the  HELP  code  was  employed  to  solve  a series  of  ! 

• impact  situations.  In  this  section  only  a brief  description  1 

of  the  numerical  techniques  and  material  models  employed  in  the 
calculations  is  presented,  since  the  HELP  code  is  well  docu- 

• mented  ( Reference;.  3 and  4)  . 

HELP  is  a two-dimensional,  multimaterial,  Eulerian  code 
for  solving  material  flow  problems  in  the  hydrodynamic  and 
elastic-plastic  regimes.  Although  the  code  is  basically 
Eulerian,  material  interfaces  and  free  surfaces  are  propagated 
through  the  c&lculational  mesh,  in  a Lagrangian  manner,  as 
discrete  interfaces  across  which  the  materials  are  not  allowed 
to  interdiffU3e,  This  interface  treatment  gives  HELP  a dis- 
tinct advantage  over  other  pure  Eulerian  codes  and  all  i it 
to  be  applied  to  the  solution  of  a variety  of  complex,  ilti- 
material  problems . The  basic  HELP  techniques  have  been  exer- 
cised over  a wide  range  of  solutions  including  hypervaloclty 
and  ballistic  impact,  fragmentation  munitions,  shaped  charge 
jet  formation  and  shaped  charge  jet  penetration.  Considerable  : 

confidence  has  been  gained  in  the  basic  numerical  model.  j 


2.2  MATERIAL  MODEL 

The  material  model  employed  in  HELP  includes  an  equation 
of  state,  a deviatoric  constitutive  relation  for  elastic  and 
plastic  deformations,  and  a failure  criterion.  These  will  be 
discussed  separately  in  the  following  paragraphs.  The  material 
constants  used  for  the  heavy  metal  projectile  and  for  the  steel 
target  are  given  in  Table  1. 


3.  Hageman,  L.  J.  and  J.  M.  Walsh,  "HELP,  A Multiple-Material 
Eulerian  Program  for  Compressible  Fluid  and  Elastic-Plastic 
Flows  in  Two  Space  Dimensions  and  Time,"  I and  II,  BRL-CR- 
39,  AD  Nos.  726459  and  7626460,  Ballistic  Research 
Laboratories,  Aberdeen  Proving  Ground,  Maryland,  May  1971. 

4.  Hageman,  L.  J.,  D.  E.  Wilkins,  R.  T.  Sedgwick  and  J.  L. 
Waddell,  "HELP,  A Multi-Material  Eulerian  Program  for 
Compressible  Fluid  and  Elastic-Plastic  Flows  in  Two  Space 
Dimensions  and  Time,"  Revised  Edition,  Systems,  Science 
and  Software  Report  SSS-R-75-256 4 , July  1975. 
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TABLE  1.  MATERIAL  PROPERTIES  EMPLOYED 
IN  THE  STANDARD  CALCULATION 


'''vNMATERIAL 

PROPERTY'''**. 

HEAVY  METAL  ALLOY* 

STEEL 

pq  <g/cm3) 

18.62 

7.8 

2 

A (dynes/cm  ) 

1.095  x 1012 

1.28  x 1012 

2 

B (dynes/cm  ) 

1.97  x 1012 

1.05  x 1C12 

a 

0.5 

0.5 

b 

1.53 

1.5 

a 

5.0 

5.0 

6 

5.0 

5.0 

IS  (ergs/g) 

o 

r—1 

o 

H 

X 

o 

• 

9.5  x IQ10 

Bs  (ergs/g) 

6.73  x 109 

2.44  :t  1010 

E'  (ergs/g) 

2.37  X 1010 

1.02  x 1011 

2 

G (dynes/cra  ) 

7.75  X 1011 

8.0  x 1011 

Yq  (dynes/cm2) 

1.01  x 1010 

6.0  x 109 

Y^  (dynes/cm2) 

0 

0 

2 

Y2  (dynes/cm  ) 

0 

0 

Em  <er9sAf> 

2.28  x 109 

1.3  x 1010 

S (dynes/cm2) 

2.7  x 1010 

3.8  x 1010 

Values  of  A,  G,  Y0,  and  S given  in  this  table  for 
the  heavy  metal  alloy  were  employed  in  Calcula- 
tions 1,  10,  and  18.  The  values  used  in  the 
remaining  calculations  are  given  in  Table 
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2,2.1  Equation  af  State 

The  Tillotson  (Reference  5)  equation  of  state,  modified 
to  give  a smooth  transition  between  condensed  and  expanded 
states,  was  employed.  In  Equations  (1)  through  (3),  p,  E and 
p are  pressure,  specific  internal  energy  and  mass  density, 
respectively;  n - p/pQ  - y + 1,  and  p„,  a,  b,  E , E , E',  A,  B, 
a and  8 are  constants  for  the  particular  material.  For  the 
condensed  states,  i.e.,  when  p/p  > 1,  or  for  any  cold  states, 

E < Ea  the  equation  of  state  has°the  form 


P " Pc 


Ep  + Ay  + Bp 


2 


(1) 


For  expanded  hot  states,  i.e.,  when  p/p  < 1 and  E > E',  the 
equation  of  state  has  the  form  ° 3 


p - pfi  - aEp  + 


bEp 


+ 1 


+ Aye"0 (pQ/p-1) 


e-a(po/p-l)  # (2) 


A smooth  transition  between  the  condensed  and  expanded  states 
is  insured  by  a transition  equation  for  the  intermediate  region 
defined  by  E < E < E'  and  p/pQ  < 1.  This  blended  portion  of 
the  equation  of  state  has  the  form 


P 


<E  - Es>PE  * <Es  " E)pc 
K " Es 


(3) 


2.2.2  Elastic-Plastic  Constitutive  Relation  and  Yield  Criterion 

The  deviatoric  stress  increments,  dS.  .,  are  determined 
by  using  the  elastic  relation  13 


5.  Tillotson,  J.  H. , "Metallic  Equations  cf  State  for  Hyoer- 
velocity  Impact,"  General  Atomic  Report  GA-3216,  July* 1962. 
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(4) 


dS^  ■ 2Gdc^j 


where  G is  the  modulus  of  rigidity  and  de, , are  the  increments 
of  deviatoric  strain.  When  such  an  increment  of  stress  causes 
the  von  Mises  yield  criterion, 

sijsij  - 2*2  (5) 


to  be  violated,  each  stress  component  is  proportionately  reduced 
to  bring  the  stress  normally  back  to  the  yield  surface.  A vari- 
able yield  strength 


Y - (Y0  + + Y2u2)  (1  - E/Em)  (6) 


is  defined  to  account  for  the  increase  in  strength  at  high 
pressures  and  the  decrease  of  strength  at  elevated  values  of 
the  specific  internal  energy,  E. 


2.2.3  Failure  in  Tension 

Given  the  spall  threshold,  S,  based  on  plate  slap 
experiments  and  expressed  as  a negative  pressure,  the  maximum 
distension  of  a material  is  approximated  in  the  HELP  code  by 
regarding  the  pressure  as  a linear  function  of  the  bulk  modulus, 
ignoring  higher  order  terms  in  p and  any  energy  dependence. 
Therefore,  given  the  spall  threshold,  S,  and  the  bulk  modulus 
A,  the  scalar  quantity  (p/Pg)_  associated  with  the  maximum 
distension  of  a material  fallows  from 


S - A (<P/P0)a  - l)  , (7) 

isie  f 


<p/e0>,  ■ | + 1 • (8) 


2. 2. 3.1  Dependence  of  Failure  on  Specific  Internal  Energy 

Previous  to  these  three  studies,  the  material  failure 
criterion  in  HELP  was  based  solely  on  material  distension  as 


indicated  by  the  scalar  quantity  (p/p0)&.  However,  in  the 
present  parametric  studies,  this  criterion  was  extended  to 
account  for  the  effect  of  heating  on  the  material's  ability 
to  withstand  tension.  It  was  assumed  that  the  spall  thresh- 
old approached  zero  as  the  material's  specific  internal 
energy  approached  E,.  . , the  energy  needed  to  melt  the  mate- 
rial. melt 

The  density  corresponding  to  the  melt  state,  Pmelt> 
was  determined  such  that 


melt' 


Emelt> 


If  the  material  was  between  the  cold  and  melt  states,  then 


(|+  D 


'moll 


which  is  a straight-line  interpolation  in  the  E - p plane  be- 
tween the  values  associated  with  the  coid  and  melt  sidles. 


2. 2. 3. 2 Gradual  and  Abrupt  Failure  Models 


When  the  distension,  p/p  , of  the  material  in  a computa- 
tional cell  was  less  than  (p/p  °)  , the  material  was  presumed 
to  have  failed.  In  all  but  one  of  the  calculations  in  these 
studies,  the  material  was  failed  gradually  over  a characteristic 
time.  The  characteristic  time  was  determined  for  each  cell  by 
dividing  the  cell  diagonal  by  the  sound  speed  of  the  material. 
The  material  in  a cell  had  to  be  subjected  to  conditions  that 
would  cause  failure  for  a time  duration  equal  to  the  charac- 
teristic time  in  order  for  the  material  in  that  cell  to  be  com- 
pletely failed.  If  the  stress  conditions  changed  before  the 
cell  completely  failed,  the  material  was  considered  only  par- 
tially damaged.  In  order  to  keep  track  of  this  damage,  a 
failure  array,  F.  , was  added  to  the  code.  When  Fk  ■ 0 , the 
material  in  the  Kth  cell  had  not  failed,  and  whenKF.  • 1,  it 
had  completely  failed.  Intermediate  values  of  F would  be  asso- 
ciated with  cells  which  had  satisfied  the  failure  criterion  for 
a time  duration  less  than  the  characteristic  time.  The  total 
stresses  calculated  in  each  damaged  cell  were  reduced  by  the 
factor  (1  - F). 
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In  the  abrupt  failure  model  employed  in  Calculation  10, 
the  valuo  of  F was  set  immediately  to  1 once  the  material  was 
subjected  to  conditions  which  violated  the  failure  criterion. 

In  both  tho  gradual  and  abrupt  failure  models,  the 
failure  parameter,  F,  was  transported  so  as  to  associate  it 
with  a volume  of  material  rather  than  with  a fixed  computa- 
tional cell.  Therefore,  values  of  F between  0 and  1 could 
result  from  material  that  was  gradually  failing  over  a char- 
acteristic time  when  the  gradual  failure  model  was  employed, 
or  from  failed  material  being  transported  when  either  model 
was  employed. 
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SECTION  III 


DESCRIPTION  OF  THE  PARAMETRIC  STUDIES 


The  HELP  calculations  performed  in  the  two  previous 
investigations  and  in  the  current  parametric  study  involved 
the  normal  penetration  of  a heavy  metal  projectile  into  5.08- 
cm-thick  rolled  homogeneous  armor  (RHA)  at  a velocity  of 
5.5  x 101*  cm/sec.  The  material  properties  and  design  of 
the  penetrator  in  Calculation  1 of  the  first  study  was  chosen 
as  the  standard  for  the  subsequent  material  property  para- 
metric studiai  . The  material  properties  of  the  projectile  in 
Calculation  1 were  chosen  as  being  representative  of  a real 
heavy  metal  alloy.  The  dimensions  of  the  projectile  are 
given  in  Figure  i.  The  material  properties  of  the  stoel  tar- 
get were  held  constant  throughout  the  studies.  (The  material 
constants  usee,  for  tie  standard  projectile  and  the  steel  tar- 
get are  given  ir.  Table  1.) 

In  this  report  the  results  of  the  calculations  from  the 
second  study  were  combined  with  those  of  the  current  investi- 
gation so  that  firmer  conclusions  regarding  the  effects  of 
penetrator  material  proparties  on  the  penetration  process 
could  be  made.  Table  2 provides  the  calculational  matrix. 

Meat  of  the  calculations  shown  in  Table  2 involved  variations 
from  the  standard  valuas  assigned  to  the  spall  threshold, 
shear  yield  strength,  bulk  modulus,  and  rigidity  modulus  in 
Calculation  1.  Two  of  the  calculations,  numbers  10  and  18, 
were  performed  to  assess  the  effects  of  specific  numerical 
tachniguesi  abrupt  failure  and  tha  addition  of  an  artificial 
viscoaity  tsrm  to  cell  boundary  atrassas.  An  analysis  of  the 
results  of  Calculations  10  and  18  is  given  in  Section  IV  of 
this  report.  It  should  be  noted  here,  however,  that  the  re- 
sults of  those  calculations  indicated  that  Calculations  1 and 
2 should  be  omitted  from  the  material  property  parametric 
study  and  that  Calculation  18  should  becoma  the  standard  in 
place  of  Calculation  1. 

Table  3 indicates  the  specific  calculations  which  were 
chosen  for  investigating  the  effects  of  varying  each  of  the 
four  material  properties.  Thesa  affects  are  discussed  in 
detail  in  Section  V. 
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Figure  1.  Dimensions  of  the  30  m Heavy  Metal  Penetrator 
Employed  in  the  Material  Parametric  Studies 


Artificial  Vtico*ity  Addcii  to  t>! I Poaiw'iry 
Abrupt  Failure  Aidtl  Frployed 

Constant  and  Equal  fco  That  r»f  Ctlcik«Uon  I* 


TABLE  3.  SPECIFIC  CALI 

each  op  urn  i 

MATERIAL  PROPERTY 
Shear  Yield,  Y 

o 

Spall  Threshold,  S 
Shear  Modulus,  G 
Bulk  Modulus,  A 


ILATIONS  USED  TO  INVESTIGATE 
IUR  MATE1UAL  PROPERTY  EFFECTS 


i 


CALCULATION  NUMBERS 


3,  4,  11,  12,  13 
6,  7,  18,  19 
8,  14,  IS,  18 
16,  17,  18 
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SECTION  IV 


EFFECTS  OF  CHANGES  TO  THE  NUMERICAL  MODEL  AND 
COMPARISON  WITH  EXPERIMENT 


In  the  first  study  (ReferJJice  1) , the  design  which  was 
chosen  for  calculation  1 was  tested  experimentally  and  the 
resulting  crater  dimensions  were  compared  to  those  predicted 
by  the  HELP  code.  The  discrepancy  of  10  percent  was  somewhat 
surprising  since  the  results  of  the  HELP  code  in  other  impact 
studies  had  been  significantly  closer  to  experimental  results 
(References  6 through  8).  In  reviewing  the  precise  model  used 
in  that  first  study  (Calculations  1 and  2) , it  wa3  concluded 
that  the  addition  of  an  artificial  viscosity  term  in  the  form  of 
a cell  boundary  stress  was  not  only  unnecessary , but  it  was  prob- 
able that  by  smearing  the  stress  field  the  added  viscosity  had 
significantly  reduced  the  interface  velocity  and  the  predicted 
final  crater  depth.  The  addition  of  an  artificial  viscosity 
term  was  based  on  studies  made  on  a very  brittle  material  in 
which  small  noises  in  the  stress  field  led  to  unrealistic  prop- 
agation of  failure.  Initially , it  was  thought  that  such  unreal- 
istic failure  might  occur  in  the  heavy  metal  penetrator  even 
though  it  was  relatively  ductile.  In  the  subsequent  studies, 
however,  the  additional  artificial  viscosity  term  was  emitted 
from  most  of  the  calculations,  and  no  evidence  of  unrealistic 
failure  propagation  was  observed. 

In  the  current  study  the  first  calculation  was  repeated; 
it  is  reported  here  as  Calculation  18.  The  crater  depth  of 
Calculation  1 (with  viscosity)  was  1.382  cm,  whereas  the 


/ 
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crater  depth  of  Calculation  18  (without  viscosity)  was 
1.53  cm,  an  increase  o£  10.7  percent  which  brings  the  calcu- 
lated results  within  a percent  of  the  experimental  results. 

Crater  depth  as  a function  of  time  is  plotted  for  Calcula- 
tions 1 and  18  in  Figure  2.  As  these  curves  indicate,  the 
crater  formation  process  is  nearly  identical  in  the  two  cal- 
culations for  the  first  5 0 ys  after  impact,  From  then  on, 
however,  the  added  artificial  viscosity  in  Calculation  1 
significantly  reduces  the  rate  of  crater  growth  and  the  * 

final  crater  depth. 

Because  of  its  close  agreement  with  experiment,  Calcu-  < 

lation  18  becomes,  in  this  current  study,  the  basic  calcula- 
tion in  place  of  Calculation  1.  Calculation  2 is  omitted 
from  the  comparisons  in  this  report  since  it  also  was  part  of 
the  first  study  and  reflects  the  use  of  an  additional  artifi- 
cial viscosity  term.  Furthermore,  the  rigidity  modulus  used 
in  Calculation  2 was  10  percent  lower  than  that  used  in  Cal- 
culations 6 and  7,  thereby  making  the  previous  conclusions  on 
the  effect  of  varying  the  shear  yield  strength  unnecessarily 
tentative. 

In  the  second  study  an  additional  artificial  viscosity 
tsrm  was  employed  in  Calculation  9 in  which  the  projectile 
material  bulk  modulus  was  twice  thi  standard  value.  In  view 
of  the  above  observations  on  the  effects  of  this  added  vis- 
cosity, Calculation  9 is  omitted  from  the  current  analysis  of 
the  effects  of  varying  the  projectile  material  bulk  modulus. 
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Figure  2.  Crater  Depth  Versus  Time  for 
Calculations  1 and  18 
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SECTION  V 

RESULTS  OF  HELP  CALCULATIONS 


5.1  SUMMARY  PLOTS  AND  TABLES 

The  HELP  calculations  reported  here  followed  the  pene- 
tration process  until  the  depth  of  the  crater  became  constant. 
The  effectiveness  of  each  projectile  material  modeled  is  mea- 
sured primarily  in  terms  of  the  final  crater  depth.  The  values 
of  crater  radius  and  projectile  length  at  that  final  time  as 
well  as  the  total  energy  coupled  to  the  target  are  presented 
as  indicators  of  differences  in  the  penetration  process  and  are 
not  meant  to  represent  the  ultimate  values  of  these  parameters 
which  may  continue  to  change  after  the  final  crater  depth  has 
been  established. 

To  illustrate  the  intermediate  stages  of  the  penetration 
process/  plots  of  the  projectile/targat  configurations  at 
various  times  for  the  standard  calculation  (number  18  in  Table 
2)  are  shown  in  Figure  3.  (The  configuration  plots  for  all 
nineteen  HELP  calculations  as  well  as  the  grid  zoning  employed 
in  all  of  the  calculations  are  presented  in  Appendix  A.) 

The  results  of  the  nineteen  calculations  are  summarized 
in  Tables  4/  5,  and  6.  Table  4 lists  the  quantitative  predic- 
tions regarding  the  final  crater  depth  and  the  corresponding 
crater  radius  and  projectile  length.  Table  5 normalizes  these 
quantities  to  Calculation  18,  the  basic  case  in  the  parametric 
study.  Table  6 Indicates  the  final  time  for  each  calculation 
(i.e.,  the  time  at  which  the  final  crater  depth  was  estab- 
lished) and  the  total  energy  coupled  to  the  target  from  the 
time  of  impact  to  the  final  time. 


5.2  EFFECTS  OF  VARYING  THE  PROJECTILE  SHEAR  STRENGTH 

The  effects  of  increasing  the  shear  yield  strength  of 
the  projectile,  as  predicted  by  the  HELP  calculations/  can  be 
seen  by  comparing  the  results  of  Calculations  12,  11,  13,  3 
and  4 in  which  the  projectile  yield  strength  was  0.75,  1.51, 
2.55,  7.56,  and  15.1  Jcilobars,  respectively.  The  final  con- 
figurations of  these  calculations  are  given  in  Figure  4.  It 
is  clear  from  these  configuration  plots  that  the  hardness  of 
the  projectile  material  significantly  affects  the  final  depth 
and  shape  of  the  crater.  From  the  curves  drawn  through  the 
data  points  in  Figure  5,  it  is  apparent  that,  as  the  shear 
yield  strength  of  the  projectile  is  increased,  the  crater 
radius  decreases  almost  linearly,  whereas  the  crater  depth 
is  essentially  unaffected  by  increasing  the  projectile  shear 
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t ■ 60  ys 


t ■ 160  us 


Figure  3.  Projectile/Target'  Configuration  at  Various  Times  for 
Calculation  18  (S  - 2.7  x 1018,  Y « 1,01  x 1010,  G - 7.75  x 10 1 1 , 
A - 1.095  x 1012)  ° 
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TABLE  6.  TIME  OF  CRATER  ARREST  AND  TOTAL  ENERGY  COUPLED 
TO  TARGET  FOR  CALCULATIONS  1 THROUGH  19 


CALCULATION 

NUMBER 

TIME  OF 
CRATER  ARREST 
(ys) 

TOTAL  ENERGY 
COUPLED  TO  TARGET 
(ergs  x 10  “) 

1 

150 

1.486 

2 

160 

1.261 

3 

160 

1.415 

4 

200 

1.619 

5 

210 

1.943 

6 

235 

2.607 

7 

220 

1.764 

8 

200 

1.637 

9 

150 

1.356 

10 

150 

1.315 

11 

200 

1.793 

12 

205 

1.930 

13 

190 

1.645 

14 

210 

1.536 

15 

160 

1.629 

16 

210 

1.516 

17 

220 

1.567 

18 

160 

1.502 

19 

160 

1.544 

20 
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Figure  5a.  Plot  of  Final  Crater  Radius  Versus  Projectile 
Shear  Yield  Strength  for  Calculations  3,  4,  11,  12  and  13 
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SHEAR  YIELD  STRENGTH  OF  PROJECTILE  (kb) 


Figure  5b.  Plot  of  Final  Crater  Depth  Versus  Projectile 
Shear  Yield  Strength  for  Calculations  3,  4,  11,  12  and  13 
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yield  strength  from  0.75  to  7.56  kilobnrs  and  is  increased 
by  24  percent  when  the  projectile  shear  yield  strength  is 
increased  from  7. 56  to  15.1  kilobars. 


The  second  parameter ic  investigation  of  heavy  metal 
penetrators  predicted  a minimum  in  the  crater  depth  versus 
projectile  shear  yield  strength  curve  corresponding  to  the 
results  of  Calculation  2.  The  additional  calculations  com** 
plated  as  part  of  the  present  study  (numbers  12  and  13)  in- 
dicate a slight  upturn  in  the  curve  at  very  low  values  of 
the  projectile  shear  yield  strength.  However,  by  excluding 
the  results  of  Calculation  2 from  the  present  analysis,  the 
importance  end  magnitude  of  the  minimum  point  are  greatly 
diminished. 

Additional  data  extracted  from  these  HELP  calculations 
suggest  e rather  abrupt  change  in  the  extent  of  projectile 
deformation  aa  tha  projectile  yield  strength  is  raised  above 
7 kilobars.  it  is  already  clear  from  the  configuration  plots 
in  Figure  4 that  the  hardness  of  the  projectile  significantly 
affects  how  much  of  the  projectile  mesa  becomes  pert  of  the 
ejects  and  how  much  remains  essentially  undeformed.  In  Figure 
6 the  curve  showing  tha  relationship  between  projectile  hard- 
ness and  final  projectile  length  drops  off  abruptly  between 
yield  strengths  of  7.56  and  2.55  kilobars.  Outside  of  this 
range,  however,  the  elope  of  the  curve  predicts  a much  weaker 
dependence  of  final  projectile  length  on  projectile  hardness. 
There  appears  to  be  a relatively  narrow  range  of  values  for 
projectile  shear  yield  strength  (between  2.55  and  7.56  kilo- 
bare  in  tha  impact  situation  modeled  by  these  calculations) 
which  defines  e transition  between  the  projectile  materials 
in  which  ell  or  moet  of  the  projectile  mass  becomes  part  of 
the  ejecta  end  those  in  which  e significant  portion  of  the 
back  of  the  projectile  remains  essentially  undeformed.  This 
critical  range  of  values  probably  is  a function  of  the  impact 
velocity,  the  materiel  properties  of  the  target,  and  the  other 
materiel  properties  of  the  projectile.  For  example,  e lower 
tensile  failure  threshold  in  the  projectile  or  e higher  yield 
strength  in  the  target  might  cause  this  transition  to  occur 
at  higher  values  of  projectile  yield  strength. 

A basic  difference  in  tho  cratering  process  due  to  the 
uae  of  herd  rather  then  soft  projectile  materials  also  is  in- 
dicated by  the  crater  depth  versus  time  curves  in  Figure  7 . 

When  the  herder  projectile  materials  (Yp  > 7.56  kb)  are  used, 
Calculations  3 and  4 predict  that  the  final  crater  depth  is 
established  earlier  and  more  abruptly  than  when  the  softer 
materials  (Yg  <2.55  kb)  are  chosen,  ae  in  Calculations  11, 

12,  end  13.  By  taking  the  first  derivative  of  these  curves,  the 
velocity  of  the  projectile/target  interface  can  be  obtained. 
Figure  7 therefore  indicates  that  in  Calculations  3 and  4 the 
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FINAL  PROJECTILE  LENGTH 


1.75 


Figure  7.  Crater  Depth  Versus  Time  for 
Calculations  3,  4,  11,  12,  and  13 


interface  velocity  remains  nearly  constant  until  about  100  ps 
and  then  over  the  next  40  us  quickly  approaches  zero.  By 
contrast,  in  Calculations  11  and  12,  where  Y < 1.51  kb,  the 
interface  velocity  drops  off  very  slowly  for  175  ps  and  goes 
to  zero  only  when  the  residual  length  of  the  projectile  ap- 
proaches zero  (as  seen  in  Figure  4),  Calculation  13,  in  which 
Yg  ■ 2.55  kbar,  represents  a transition  between  the  extremes. 
Tne  deceleration  of  the  projectile/target  interface  occurs  more 
rapidly,  as  it  does  when  the  hard  materials  are  employed,  but 
the  interface  continues  to  move  until  most  of  the  projectile 
material  has  been  consumed,  which  is  characteristic  of  the 
softer  materials. 

This  difference  in  crater  growth  rate  and  in  reiiidual 
projectile  length  suggest  that  whan  the  yield  strength  is 
sufficiently  high  (^  7.56  kb),  the  projectile  is  stopped  by 
stresses  within  the"“projectile  as  well  as  by  those  within 
the  target.  Furthermore,  when  the  yield  strength  of  the  pro- 
jectile material  is  sufficiently  low  (<  1.51  kb) , these 
stresses  are  not  sufficient  to  stop  the  projectile  until  it 
is  entirely  consumed  by  the  cratering  process. 

A difference  in  the  effect  of  the  hard  and  soft  pro- 
jectile materials  is  reflected  elso  by  Figure  8 in  which  the 
total  energy  coupled  to  the  target  is  plotted  versus  projec- 
tile shear  yield  strength.  The  softer  projectiles  (Y0  £ 2.55 
kb) , which  form  wider  craters  and  thereby  distribute  their 
load  over  a larger  surface  area,  couple  more  total  energy  to 
the  target  than  do  the  harder  projectiles  (Y  > 7.56  kb). 

These  results  are  not  surprising  in  light  of  tKe  fact  that 
more  of  the  hard  projectile's  kinetic  energy  is  converted  to 
heat  by  plaatic  work  and  therefore  lees  of  it  is  available  to 
be  coupled  to  the  target.  In  addition,  the  hard  projectile, 
because  it  flows  leas  radially,  focuses  its  load  over  a 
smeller  surface  area  and  thereby  forme  a deeper,  narrow 
crater  even  though  it  couples  leas  total  energy  tc  the  target. 
The  crater  profiles  in  Figure  9 show  the  change  in  crater 
shape  as  the  yield  strength  of  the  projectile  material  is 
increased . 

To  summarize,  the  results  from  Calculations  3 and  4 
predict  e 24  percent  gain  in  crater  depth  whan  the  projectile 
materiel  shear  yield  strength  is  increased  by  a factor  of  two 
from  7.56  to  15.1  kilobars.  However,  when  this  parameter  is 
increased  from  0.75  to  7.56  kilobars,  a factor  of  ten.  Calcu- 
lations 12  end  3 predict  narrower  but  not  deeper  craters. 

Crater  depth,  therefore,  has  a nonlinear  dependence  on  projec- 
tilo  hardness  over  the  range  of  values  considered  here  (0.75 
to  15.1  kb).  Furthermore,  above  a critical  value  of  projec- 
tile yield  strength,  the  cratering  process  apparently  changes 
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TOTAL  ENERGY  COUPLED 


Figure  8.  Total  Energy  Coupled  to  the  Target  Versus 
the  Shear  Yield  Strength  of  the  Projectile  for 
Calculations  3.  4,  11.  12.  and  13 
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character;  i.e.,  the  residual  projectile  length  . ceases 
rapidly,  the  growth  of  the  crater  end/,  abruptly,  and  before 
all  of  the  projectile  is  consumed.  This  critical  value  is 
very  probably  a function  of  the  impact  velocity,  the  target 
material  properties,  the  other  projectile  material  proper- 
ties, and  the  initial  dimensions  of  .e  projectile. 


5.3  EFFECTS  OF  VARYING  THE  TENSILE  FAILURE  THRESHOLD  OF 

THE  PROJECTILE 

The  effects  of  changing  the  tensile  failure  threshold 
of  the  projectile  can  be  seen  by  comparing  the  results  of 
Calculations  6 and  7 and  those  of  Calculations  IS  and  19. 

As  indicated  in  Table  2,  the  projectiles  in  Calculations  6 
and  7 were  relatively  soft,  having  a yield  strength  of  3.78 
kilobars,  whereas  the  projectiles  in  Calculations  18  and  19 
had  a yield  strength  of  10.1  kilobars,  representing  a rela- 
tively hard  material.  In  both  aata  the  spall  threshold  of 
one  was  27  -kilobars  (Calculations  7 and  18) , and  the  other 
was  18  kilobars  (Calculations  6 and  19) . 

The  failure  criterion  employed  in  these  calculations 
is  described  in  Section  XI.  However,  a brief  description  of 
how  the  criterion  is  applied  seems  appropriate  at  this  point. 
When  the  material  in  a calculational  cell  satisfies  the 
failure  criterion  over  a characteristic  time,  the  code  presumes 
that  the  material  has  become  a collection  of  small  particles 
which  is  unable  to  withstand  tension  (negative  pressure)  and 
has  no  strength  in  shear.  However,  the  code  does  assume  that 
this  failed  material  can  later  sustain  hydrostatic  pressures 
if  it  becomes  sufficiently  compressed  or  heated. 

As  the  configuration  plots  and  the  crater  depth  versus 
time  plots  in  Figures  10  through  12  indicate,  varying  the 
tensile  threshold  of  a soft  projectile  (YQ  ■ 3.78  kb)  has  a 
significant  effect  on  the  final  crater  geometry;  whereas, 
this  parameter  when  varied  from  18  to  27  kilobars  in  a hard 
penetrator  material  (YQ  - 10.1  kb)  has  essentially  no  effect 
on  the  final  crater  dimensions  or  on  the  penetration  process. 
These  conclusions  are  summarized  by  Figures  13a,  13b,  and  13c 
in  which  final  crater  depth,  crater  radius,  and  projectile 
length  are  plotted  versus  the  projectile  spall  threshold  for 
Calculations  6 and  7 and  Calculations  18  and  19.  For  the 
hard  projectiles  (18  and  19)  all  of  these  quantities  are 
essentially  the  same.  For  the  soft  projectiles  (6  and  7), 
however,  crater  depth  is  clearly  affected  by  changes  in  the 
tensile  failure  of  the  projectile;  the  crater  radius  and  the 
final  length  of  the  projectile  are  also  affected  but  to  a 
lesser  degree.  In  fact,  the  final  length  of  the  projectile 
is  more  sensitive  to  projectile  yield  strength  than  it  is  to 
spall  threshold  as  is  evidenced  by  the  spread  in  the  data  for 
the  soft  and  hard  projectiles  in  Figure  13c. 
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Figure  11.  Crater  Depth  Versus  Time  for  Calculations 
6 and  7 (Soft  Projectiles) 
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Figure  13.  Final  Crater  Depth,  Crater  Radius,  and 
Length  of  Projectile  Versus  Spall  Threshold  for 
Calculations  6,  7,  13,  and  19 
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These  results  raise  two  interesting  questions.  (1) 

Why  does  decreasing  the  spall  threshold  of  a soft  projectile 
make  it  a better  penetrator?  (2)  Why  doee  decreasing  the 
spall  threshold  of  a hard  projectile  have  no  effect  on  its 
penetrability?  It  appears  that  the  answers  to  both  questions 
involve  the  effect  of  projectile  shear  yield  strength  on  the 
penetration  process. 

In  the  first  instance*  the  increased  tensile  failure  in 
the  projectile  has  two  effects  or.  the  cratering  process.  It 
drastically  reduces  the  plastic  work  of  the  projectile  mate- 
rial* and  it  narrows  down  the  area  across  which  the  projectile 
load  is  delivered.  The  first  of  those  effocts  enables  the  pro- 
jectile to  couple  more  of  its  energy  to  the  target*  and  the 
second  enables  it  to  form  a deeper  crater.  Figures  14a  and 
14b  show  the  effect  of  projectile  spall  threshold  on  the  plas- 
tic work  of  the  projectile  and  the  total  energy  coupled  to  the 
target. 

In  the  second  instance,  it  can  be  concluded  that  the 
high  strength  of  the  hard  projectile  inhibits  the  radial 
flow  of  the  material  within  the  crater  volume  which*  in  turn, 
prevents  hydrostatic  tensions  above  18  kilobars  from  occur- 
ring. This  conclusion  is  supported  by  Figure  15  which  com- 
pares the  extent  of  failed  material  at  160  ys  for  the  cal- 
culations involving  soft  and  hard  projectiles,  both  having 
a spall  threshold  of  18  kilobars.  In  Figure  15  the  shaded 
cells  contain  material  which  is  expanded  (p  < pQ)  and  which 
at  some  point  has  failed.  The  calculations  proaict  consider- 
able failure  of  projectile  material  within  the  crater  volume 
for  the  soft  projectile  (Calculation  6)  and  essentially  none 
for  the  hard  projectile  (Calculation  19).  It  can  be  concluded* 
therefore*  that  the  spall  threshold  for  the  hard  projectile 
must  be  lower  than  18  kilobars  for  it  to  affect  the  projec- 
tile's penetrability. 

By  comparing  the  crater  depth  curves  versus  time  for 
Calculations  6 and  7 (Figure  11)  and  for  Calculations  3*  4* 

11*  12,  and  13  (Figure  7),  it  is  apparent  that  the  yield 
strength  of  the  projectile  does  affect  the  early  stages  of 
the  cratering  process*  whereas  the  ductility  of  the  projec- 
tile does  not  have  a significant  effect  on  the  crater  depth 
until  after  100  ys.  These  results  seem  to  suggest  that  to 
optimize  the  performance  of  a heavy  metal  penetrator  against 
a thick  steel  target  it  is  important  for  it  to  have  a high 
shear  yield  strength  for  the  early  stages  of  penetration  and 
to  have  a low  spall  threshold  for  the  late  stages. 


5.4  EFFECTS  OF  VARYING  THE  SHEAR  MODULUS  OF  THE  PROJECTILE 

By  comparing  the  results  of  Calculations  3*  14*  15, and 
18,  the  effects  of  varying  the  projectile  material  shear 
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Figure  15.  Configuration  Plot#  of  Calculations  6 and  19 
at  160  ys.  Failed  Cells  With  Zero  Stresses  are  Shaded. 
The  Spall  Threshold  in  Both  Calculations  was  18  Kilobars 
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modulus  can  be  observed.  The  ahaar  modulus  of  a real  heavy 
natal  alloy,  ao  modeled  by  Calculation  18,  is  7.75  x 1C11 
dynea/cra2.  This  value  halved,  i.e.,  3.875  x 10* 1 dynes/cm2, 
was  used  in  Calculation  15,  and  in  Calculations  8 and  14  the 
standard  value  was  increased  by  factors  of  two  and  four  to 
give  shear  moduli  of  15.5  x 101 1 dynes/cm2  and  31.0  x 10l 1 
dynes/cm2,  respectively , Those  four  calculations  therefore 
offer  results  over  a substantial  range  of  projectile  shear 
moduli. 

As  the  configuration  plots  in  Figure  16  suggest,  the 
final  crater  dimensions  are  rather  insensitive  to  changes  in 
the  rigidity  of  the  projectile.  All  of  the  craters  are  nearly 
hemispherical,  and  there  is  less  than  a 6 percent  variation  in 
the  final  crater  depth  or  radius  from  the  average  values 
(15  ■ 1.617  cm,  5 ■ 1.514  cm}  predicted  by  these  four  calcula- 
tions. 

The  final  crater  depths  and  radii  ari  plotted  versus 
projectile  shear  modulus  in  Figure  17a  and  17b.  Aside  from 
the  fact  that  there  is  some  scatter  in  the  crater  depth  data, 
a nearly  horizontal  line  can  be  drawn  through  both  the  crater 
depth  and  crater  radius  data.  Likewise,  as  indicated  by 
rigure  17c,  essentially  the  same  total  energy  is  coupled  to 
the  target  during  the  cratering  process  by  all  four  projec- 
tiles. Therefore,  these  results  clearly  predict  a very  weak 
relationship  between  projectile  rigidity  and  penetrability. 

The  final  length  of  the  projectile,  on  the  other  hand, 
appwars  to  increase  significantly  as  the  projectile  shear 
modulus  drops  below  15.5  x 10u  dynes/cm1,  as  seen  in  Figure 
18a.  It  should  be  noted  again,  however,  that  the  final  pro- 
jectile length  is  obtained  at  the  time  in  the  calculation 
that  the  final  crater  depth  is  established.  As  a result, 
these  data  points  arc  an  indication  of  how  much  of  the  pro* 
jectile  material  is  consumed  in  establishing  a final  crater 
depth,  and  they  are  not  necessarily  an  indication  of  the  pro- 
jectile length  when  all  of  the  projectile's  kinetic  energy 
has  been  dissipated.  These  calculations  seem  to  indicate 
that  whila  the  back  of  the  projectile  is  still  moving  toward 
the  target,  a significant  amount  of  the  projectile  material 
flows  radially  and  becomes  part  of  ^he  lip  and  debris  before 
coming  close  enough  to  the  projectile/target  interface  to 
affect  the  depth  of  the  crater.  In  the  plot  of  residual  pro- 
jectile kinetic  energy  versus  projectile  shear  modulus  in 
Figure  18b,  it  is  evident  that  the  projectiles  with  the  low- 
est shear  modulus  have  the  largest  residual  kinetic  energy 
as  well  as  the  greatest  length  when  the  crater  depth  is 
finalised.  Consequently,  at  the  time  that  the  projectile 
kinetic  energy  goes  to  zero  in  Calculations  15  and  18,  the 
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Figure  17.  Crater  Depth.  Crater  Radius,  and  Total 
Energy  Coupled  to  the  Target  Versus  Time  for 
Calculations  8.  14.  15.  and  16 
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projectile  length  might  be  reduced  to  that  of  the  projectiles 
in  Calculations  8 and  14.  Therefore , it  is  possible  that  the 
true  final  dimensions  of  the  projectile  as  well  as  those  of 
the  crater  are  essentially  unaffected  by  the  rigidity  of  the 
projectile. 

Finally,  the  similar  shape  and  slope  of  the  crater 
depth  versus  time  curves  in  Figure  19  for  Calculations  8,  14, 

15,  and  18  suggest  that  the  growth  rate  of  the  crater  is  also 
essentially  independent  of  the  rigidity  modulus  of  the  pro- 
jectile material. 

5.5  EFFECTS  OF  VARYING  THE  BULK  MODULUS  OF  THE  PROJECTILE 

Calculations  16,  17,  and  18  were  performed  for  the  pur- 
pose of  investigating  the  effectfe^f  varying  the  bulk  modulus 
of  the  projectile  material.  In  the  standard  calculation  (18) 
the  bulk  modulus  of  a real  heavy  metal  alloy,  1.095  x 101Z 
dynes/cm2,  was  used.  One-half  of  this  value,  or  5. 975  x 10* 1 
dynes/cm1,  was  used  in  Calculation  16  and  four  times  the 
standard  value,  or  4.38  x 1012  dynes/cm2,  was  used  in  Calcu- 
lation 17. 

Figure  20  shows  the  predicted  final  projectile/target 
configurations  for  these  three  impact  situations.  It  is 
evident  from  the  projectile  length  and  lip  formation  at 
crater  arrest  that  the  cratering  process  is  somewhat  affected 
by  the  bulk  modulus  of  the  projectile.  However,  the  plots  in 
Figure  21  show  that  the  final  dimensions  of  the  crater  and 
the  energy  coupled  to  the  target  are  relatively  Insensitive 
to  the  bulk  modulus  of  the  projectile  within  the  range  repre- 
sented here.  It  is  apparent  from  the  arater  depth  versus 
time  curves  in  Figure  22  that  the  growth  rates  of  the  craters 
in  these  three  calculations  are  not  strongly  affected  by 
changes  in  the  projectile  bulk  modulus.  The  curves  remain 
coincident  or  parallel  until  75  us,  and  they  all  begin  to 
level  off  at  about  125  u*»  Only  the  very  late  stages  of  the 
penetration,  after  150  us,  do  their  slopes  significantly 
diverge . 

As  indicated  in  Table  4,  the  spall  threshold,  S,  of  the 
projectile  in  Calculations  16  and  17  was  different  from  the  27- 
kilobar  threshold  used  in  Calculation  18.  This  difference  arose 
from  using  the  same  value  of  (p/p  ) for  the  projectile  in  the 
three  calculations  despite  the  change  in  the  material  bulk 
modulus,  A,  (Note  from  subsection  2.2.3,  (p/p  )-  ■ M 1). 

The  analysis  of  the  effect  of  varying  the  projeotilaspall 
threshold  in  subsection  5.3  reveals,  however,  that  when  hard 
material  projectiles  are  used,  the  cratering  process  and  the 
final  crator  dimensions  are  essentially  unaffected  by  a 50 
percent  increase  in  the  spall  threshold  (from  18  to  27  kb) . 


TIME  (US) 


Figure  19.  Crater  Depth  Versus  Time  for  Calculations 
8,  14,  15,  and  18 
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Figure  20.  Predicted  Final  Configurations  of 
Calculations  10 , 17,  and  18 


I 


TOTAL  ENERGY  COUPLED 
TO  TARGET  (ergs) 


0 25  50  75  100  125  150  175  200  225 

TIME  (us) 


Figure  22.  Crater  Depth  Versus  Time  for 
Calculations  16,  17,  and  18 
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SECTION  VI 


CONCLUSIONS 


In  the  current  study  the  HELP  code  was  employed  to 
calculate  eight  impact  situations  to  augment  two  earlier  para- 
metric investigations  involving  the  effects  of  projectile 
material  properties  on  penetration  into  a thick  steel  target. 
In  all  of  the  calculations  heavy  metal  penetrators  were  im- 
pacted normally  into  a 5 . 08-cm-thick  steel  plate  at  a velocity 
of  5.5  x 10• **  cm/sec.  The  projectile  material  properties 
that  were  varied  were  the  shear  yield  strength,  the  tensile 
failure  threshold,  the  shear  modulus,  and  the  bulk  modulus. 
Analysis  of  the  results  from  the  current  calculational  effort 
combined  with  the  results  from  the  two  previous  studian  has 
led  to  the  following  conclusions-  These  conclusions  apply  to 
the  particular  impact  situation  studied. 


• There  exists  a critical  value  of  projectile 
material  shear  yield  strength  above  which  the 
crater  depth  is  significantly  increased  by 
hardening  the  projectile  material.  Below  this 
critical  value,  the  crater  depth  is  not  strongly 
affected  by  increasing  the  shear  yield  strength 
of  the  projectile,  and  the  trend  is  actually  re- 
versed when  very  soft  (Y0  < 1.0  kb)  projectile 
materials  are  employed. 

• Within  the  range  of  projectile  material  shear 
yield  strengths  considered  in  these  studies 
(0.75  kb  < Yg  1 15.1  kb) , the  softer  projec- 
tiles (yq  <.  2.55  kb)  form  wider  craters  and 
couple  more  total  energy  to  the  target  than  do 
the  harder  projectiles  (7,56  kb  < Y0  < 15.1  kb). 

• Since  the  crater  depth  resulting  from  an  impact 
of  a soft  projectile  (Yc  < 2.55  kb)  is  not 
established  until  nearly  all  of  the  projectile 
mass  has  been  consumed  by  the  penetration  pro- 
cess, the  initial  dimensions  of  the  projectila 
becomes  an  important  parameter  when  the  projec- 
tile material  is  relatively  soft. 


• Increasing  the  brittleness  of  a low  yield  strength 

projectile  material  (Yc  « 3.78  kb)  by  decreasing 
the  tensile  failure  threshold  from  27  to  18  xilo- 
bara  results  in  a significantly  deeper  crater. 
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• Increasing  the  brittleness  of  a high  yield 
strength  projectile  material  (YQ  - 10.1  kb) 

by  lowering  the  tensile  failure  threshold  from 
27  to  18  kilobars  has  essentially  no  effect  on 
the  crater  dimensions.  The  10,1-kilobar  yield 
strength  of  the  projectile  material  apparently 
prevents  it  tror.  developing  tensions  that  ex- 
ceed the  lower  l:-kilobar  failure  threshold. 

• The  crater  dimensions  remain  essentially  the 
same  when  the  shear  modulus  of  the  projectile 
material  is  increased  an  order  of  magnitude 
from  3.875  x 10M  dynes/cm2  to  3.1  x 1012 
dynea/cm2 . 

• Changes  in  the  bulk  modulus  of  the  projectile 
material  between  5.975  x 1011  dynes/cm1  and 
4.38  x 10 12  dynes/cm2  do  not  strongly  affect 
the  final  crater  dimensions. 


The  results  of  these  parametric  studies  indicate  that 
the  yield  strength  and  the  tensile  failure  threshold  of  a 
heavy  metal  penetrator  strongly  influence  its  ability  to 
penetrate  a thick  steel  target.  However,  a closer  examina- 
tion of  the  interdependence  of  the  tensile  failure  threshold 
and  the  shear  yield  strength  of  the  projectile  material  is 
required  in  order  to  optimize  the  effects  of  both  character- 
istics. In  particular,  by  employing  numerical  tools  such  as 
the  HELP  code,  both  materia?  parameters  can  be  varied  in  a 
series  of  calculations  to  find  the  optimum  strength  charac- 
teristics for  a heavy  metal  penetrator. 

Perhaps  the  most  important  extension  of  these  para- 
metric studies,  however,  should  be  the  investigation  of 
oblique  impacts.  Some  of  the  above  conclusions  may  very 
well  be  altered  when  the  targets  are  impacted  with  signify 
cant  degrees  of  obliquity.  For  example,  brittleness  of  the 
projectile  material,  which  seems  to  be  an  advantage  in  the 
normal  impact  situation  studied  here,  may  actually  cause  a 
long  heavy  metal  penetrator  to  prematurely  break  in  half 
when  it  impacts  a target  obliquely.  Newly  developed  three- 
dimensional  numerical  techniques  are  available  for  the  study 
of  oblique  impacts. 

It  is  felt  that  finite  difference  numerical  methods, 
when  coordinated  with  experimental  programs  and  when  applied 
to  parametric  investigations  such  as  the  one  reported  here, 
can  efficiently  and  economically  provide  valuable  guidance 
to  the  designers  of  heavy  metal  penetrators. 
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APPENDIX  A 


PROiTECTI LE/TARGET  CONFIGURATIONS  AT  VARIOUS  TIMES 
TOR  CALCULATIONS  1 THROUGH  19 


In  Figure  A-l,  the  grid  zoning  employed  in  all  of  the 
HELP  calculations  is  shown  relative  to  the  dimensions  of  the 
target  and  projectile.  The  remaining  figures  of  this  appen- 
dix give  the  predicted  projectile/target  configurations  at 
various  times  for  the  nineteen  HELP  calculations  used  in  the 
current  analysis.  In  all  cases  calculational  times  of  5,  30, 
60,  90,  and  120  v>s  and  the  final  time  corresponding  to  ar- 
rested crater  growth  are  given.  In  each  figure  caption  the 
spall  threshold,  S,  shear  yield  strength,  Y , bulk  modulus  A, 
and  shear  modulus,  G are  indicated  in  units  of  dynes/cm*. 
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Figure  A-l.  Grid  Zoning  Employed  in  HELP  Calculation*. 
Only  the  Right  Half  of  This  Figure  Represents  the  HELP 
Grid.  The  Code  Uses  Cylindrical  Coordinates  to  Reflect 
the  Axial  Axis  of  Symmetry. 
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t « SO  US 


t ■ 60  us  t - 150  us 

Figure  A-2.  Projectile/Target  Configuration  at  Various 
Tinas  for  Calculation  1 (S  - 2.7  x 10s Y - l.oi  x 101', 
G - 7.75  x 10n,  A - 1.095  x 10,a) 
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t • 60  us 


t ■ 160  us 


Figure  A- 
Times  for 
G - 7.02 


Projectile/Target  Configurations  at  Various 
Calculation  2 (S  ■ 3.6  x 101*,  Y » 3.59  x 10*, 
: 1011  , A - 1.095  x 101  *) 
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t ■ 60  ys 


t ■ 200  US 


Figure  A-6.  Projectile/Target  Configurations  at  Various 
Times  for  Calculation  5 (S  - 1.2  x 10“,  Yc  - 5.04  x 10% 
G - 7.75  x 10u,  A - 1.095  x 10*  2) 


t ™ ?r  us 


t ■ 1Z0  us 


t ■ 60  us  t ■ 220  u» 

Figure  A-7.  Projectile/Target  Configurations  at  Various 
Times  for  Calculation  6 (S  - 1.8  x 10lt,  Y ■ 3.78  x 10*, 
G - 7.75  x 10",  A - 1.095  X 10")  ° 
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t « 60  us 


t ■ 220  US 


Figure  A-8.  Pro jactils/Target  Configurations  at  Various 
Tirass  for  Calculation  7 (S  ■ 2.7  x 101',  Y ■ 3.78  x 10* 
Q - 7.75  x 10n,  A - 1.095  x 105  *)  ° 
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* A 


Figure  A-14.  Projectile/Targot  Configurations  at  Various 
Times  for  Calculation  13  (S  ■ 3.6  x 101*,  Y ■ 2.55  x 10*, 
G - 7.75  X 10“,  A « 1.095  x 1011)  ° 
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t * 60  pi 


t - 160  pis 


Figure  A-16.  Projectile/Target  Configurations  at  Various 
Times  for  Calculation  15  (S  - 2.7  x 10l#,  Y ■ 1.01  x 10l°, 
G - 3.875  X 1011,  A - 1.095  x 101*)  ° 
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t - 60  (is 


t ■ 220 


Figure  A -18.  Projectile/Target  Cor.f  igurutions  at  Various 
Times  for  Calculation  17  (S  ■ 2.7  x 10**,  Y ■ 1.01  x 10* 
G - 7.75  X 10“,  A - 4.38  x I0,J)  ° 


t ■ 60  us 


t ■ 160  us 


Figure  A -19 . Projectile/Target  Configurations  at  Various 
Times  for  Calculation  18  (5  - 2.7  x 10l#,  Y ■ 1.01  x lO10, 
G - 7.75  x 10M,  A - 1.095  x 1011)  ° 


t - 60  ps  t ■ 160  pa 


Figure  A-20.  Projectile/Target  Configurations  at  Various 
Times  for  Calculation  19  (S  - 1.8  x 10* 4 , Yrt  » 1.01  x 101#, 
G - 7,75  X 1011,  A - 1.095  x 10 1 a ) ° 
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(The  reverse  of  this  pege  is  blank) 
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INITIAL  DISTRIBUTION 


Hq  U3AF (AF/SAMI ) 1 
Hq  AFSC/OLCAW  1 
AFIT/LD  1 
CIA-CRE/ADD/Pubs  2 
AUL/AUL-LSE-70-239  1 
Redstone  Scl  Info  Cen/Ooc  Sec  2 
Ballistic  Res  Labs  1 
SARPA-TS  1 
USN  NAV  ORD  LAB/Tech  Lib  2 
NAV  ORD  STN/Tech  Lib  1 
NAV  SYS  Cen/Tech  Lib  154  1 
DDC/TC  2 
Ogden  ALC/MNOP  2 
US  Atomic  Energy  Com/  Hq  Lib  1 
TAWC/TRADOCLO  1 
AFATL/DL  1 
AFATL/DLOSL  9 
AFATL/DLDG  1 
ASD/ENFEA  1 
ASO/CNYEHM  1 
AFIS/INTA  1 
Systems,  Science  and  Software  6 
AF  Materials  Lab/MAMD  1 
Georgia  Inst  of  Technology  3 
Lehigh  University  1 
University  of  Florida  2 
AFATL/DLJW  10 
AF  Off  of  Scl  Rseh/NAM  2 
AEDC/ARO,  Inc.  1 
Hq  TAC/DRA  1 
Hq  USAFE/DOQ  1 
Hq  PACAF/DOO  1 
AFWL/LR  2 
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